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Abstract 

We  present  an  asymptotically  optimal  algorithm  to  locate  all  the  axes  of  mirror  symmetry  of  a  planar  point 
set  The  algorithm  was  derived  by  reducing  the  2-D  symmetry  problem  to  linear  pattern-matching.  Optimal 
algorithms  for  finding  rotational  symmetries  and  deciding  whether  a  point  symmetry  exists  are  also  presented. 


1.  Introduction 

A  set  of  points  P  in  the  plane  has  an  axis  of  mirror  symmetry  A  when  for  every  point  p  of  P  not  lying  on  A 
there  is  another  point  p'  in  P  s.t  A  is  the  perpendicular  bisector  of  the  line  pp1.  P  has  a  rotational  symmetry 
of  a  when  rotating  P  about  its  centroid  C  by  a  is  an  identity  operation  on  P.  P  has  a  point  symmetry  (which 
must  be  at  C)  precisely  when  P  has  a  rotational  symmetry  of  w . 

This  note  presents  an  optimal  0(n  log  n)  algorithm  for  discovering  all  the  mirror  symmetries  of  an  n  point 
set  P  in  detail  and  describes  the  changes  needed  to  detect  the  rotational  symmetries  (and  hence  the  existence 
of  a  point  symmetry).  Lower  bounds  are  shown  for  each  problem.  Mirror  symmetries  are  the  objects  of 
interest  until  section  5. 

The  2-D  mirror  symmetry  problem  is  reduced  to  a  1-D  pattern-matching  problem  for  which  fast  solutions 
are  well-known.  Any  A  must  pass  through  the  centroid  C  of  P ,  so  the  points  are  first  translated  so  that  (0,0) 
corresponds  to  C .  After  expressing  the  points  in  polar  coordinates,  sort  them  increasing-distance-within- 
increasing  angle  from  a  reference  direction.  For  each  unique  angle  (at  most  n),  replace  the  set  of  points 
residing  at  that  angle  (i.e.,  >  1  point)  by  a  tuple  which  is  simply  a  list  of  their  distance  components.  Let  the 
number  of  unique  angles  be  m  (^  n).  Consider  the  result  as  a  length  2m  string  F,  the  symbols  of  which  are 
alternately  tuples  and  the  angles  between  adjacent  tuples.  The  mirror  symmetries  of  P  correspond  exactly  to 
the  length  2m  subsequences  of  FF  which  are  palindromes. 

The  palindromes  can  be  discovered  by  a  fast  one-dimensional  string  matching  algorithm,  looking  for 
occurrences  of  the  reversal  of  F  in  FF .  One  such  algorithm  is  that  of  Knuth.  Morris  and  Pratt  [1]  ("KMP”) 
which  permits  the  detection  of  all  occurrences  of  a  pattern  within  a  text  in  time  proportional  to  the  sum  of  the 
lengths  of  the  text  and  pattern. 

Section  2  of  this  report  contains  the  algorithm,  and  section  3  contains  proof  of  why  the  algorithm  works.  In 
section  4  we  demonstrate^that  it  is  impossible  to  improve  the  asymptotic  time  bound.  In  section  S  wq  present 
a  variant  of  the  basic  algorithm  which  detects  the  rotational  symmetries  and  point  symmetry,  and  a  proof  that 
these  are  also  optimal  ^ _ _ 
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2.  Algorithm 

Input :  n  point  planar  point  set  P.  ( p.  for  i  *-0(l)n-l). 

Output :  'ITic  centroid  of  P  and  the  orientations  of  each  of  its  axes  of  mirror  symmetry. 

1.  Find  the  centroid  C  of  P.  Translate  the  point  set  so  that  the  origin  coincides  with  C. 

1  Select  an  arbitrary  reference  direction  (for  convenience  we  choose  the  direction  of  pQ ).  Represent 
the  points  in  polar  coordinates  with  the  angle  component  as  measured  anticlockwise  from  the 
reference  direction,  denote  point  p.  as  (r.  ,8 . ). 

3.  Sort  the  points  by  incrcasing-distancc-within-increasing-anglc.  Delete  any  points  that  have  zero 
distance  (r,  =  0).  Let  the  number  of  different  angles  be  m  <  n . 

4.  For  each  of  the  m  different  angles,  represent  its  set  of  points  in  a  single  tuple .  which  simply  holds 
the  set  of  distances  at  that  angle.  Beginning  at  the  reference  direction,  proceed  through  the  list  of 
tuples  in  order  of  increasing  angle,  generating  the  length  2m  string  F :  at  the  current  position 
append  the  tuple  to  the  string;  move  to  the  next  tuple,  appending  the  angle  traversed  to  the  string, 
finish  when  returning  to  the  reference  direction.  The  first  element  of  this  string  is  a  tuple,  angles 
and  tuples  alternate  within  it.  Create  the  length  2m  string  R  and  the  length  4m-2  string  F'  as 
shown  below: 


R  ®  *0  *2m-l  f2m-2  ~  *1 

F7  *  f0  f0  *i  •••  f2m-3 

3.  Employ  a  string-matching  algorithm  such  as  KMP  to  locate  all  of  the  matches  of  R  in  F;  .  From 
the  definition  of  F  the  only  matches  are  possible  beginning  at  even  indices  (since  fQ  is  a  tuple). 
Denote  the  size  t  (0  <;  t  <;  m)  list  of  indices  of  F 1  at  which  a  match  can  begin  by  I j(j  =0(1)  t-1). 

6.  Compute  for  each  match  index  1^  the  orientation  of  an  axis  of  mirror  symmetry  A j .  Let  =  \.  /2, 


kj is  even  =>  A^  passes  through  tuple  fk  . 


kj  is  odd 


=>  A, bisects  angle  f.  . 

J  J 


3.  Details 

In  this  section  we  demonstrate  that  there  is  a  1: 1  mapping  between  the  axes  of  mirror  symmetry  of  the  point 
set  and  the  successful  matches  of  R  in  F ' .  We  also  derive  die  formula  for  finding  the  orientation  of  the 
corresponding  axis  from  a  match. 

A  palindrome  is  a  string  of  symbols  with  the  property  that  in  whichever  direction  it  is  read,  it  win  give  the 
same  string,  -  sq.i  =  Vl  ~  sl*0  ■  An  '  ‘palindrome  is  a  string  with  the  property  that  if  the  symbol  with 
index  i  ts  deleted  and  the  left  remaining  substring  is  appended  to  the  right  of  the  right  remaining  substring, 
then  the  new  string  is  a  palindrome,  s/+1 ...  s^  Sq  -  sM  =  sM ...  Sg  s^j  ».  sy  +1 . 


Because  the  tuples  do  not  necessarily  have  a  unique  distance  associated  with  them,  we  can  consider  the  set 
of  tuples  {f^ ,  for  /  =  0  to  m-1 }  as  laid  out  on  the  perimeter  of  a  circle  (diameter  unimportant,  centre  at  the 
centroid  of  the  point  set),  with  tuples  f2 .  and  f2/+2  separated  by  the  angle  f2/  +  1  •  Tuples  f2n,2  and  fQ  arc 
separated  by  angle  .  The  axes  of  this  representation  arc  precisely  the  axes  of  the  original  point  set. 

The  notation  used  is  that  of  the  algorithm  of  section  2 .  F  =  fQ  f, ...  ,  where  m  is  the  number  of  tuples 

and  an  element  f.  is  either  an  angle  (/is  odd),  or  a  tuple  (/is  even).  R  =  ffl f^  f^  ...  fj,  and 
F;  =  fQ  f, ...  fQ  ^ ...  f^  j  .  We  say  that  R  matches  F;  at  index  /  when  the  length  2m  substring  of  F; 
f. ...  f.  2  fQ ..  fM  matches  R.  For  all  such  matches  /must  be  even  because  fQ  is  a  tuple  and  the  tuples  all 
appear  with  even  index  in  F. 

lemma  For  any  k  (0  <  k  <  m) :  F  is  a  k  -palindrome  iff  R  matches  F 1  at  index  2k. 

Define  substrings. 


*1  = 

(Length  k  -1) 

h  ~  f2*+l“*f2m-l 

(Length  2m-2*  -1) 

*3  * 

(Length  k ) 

bl  ~ 

(Length  k  -1) 

b2  =  f2m-l  “  f2A  +1 

(Length  2m-2Jt  -1) 

b3  =  f2*  ~  f*  +1 

(Length  k ) 

When  F  is  a  it -palindrome  and  IK  in  we  know  that  ft  +1  „ft-1  is  a  palindrome.  More  specifically, 
aj  f^  aj  aj  =  bj  f„  b2  b3 .  Therefore  f^  =  f„ ,  aj  =  b2  .a^  =  b2 ,  and  a3  =  bj  using  the  substring  lengths 
given  above. 

When  R  matches  F;  at  index  2k  we  know  ^  a3  fk  a1  =  fQ  b2  bj  fk  b2 .  The  equivalences  here  are  again 

^2Jk  =  ^0  ’  ai  =  bl  ,a2  =  b2  ’  *3  =  b3  • 

Because  the  equivalences  enforced  are  the  same  in  both  eases  the  two  situations  are  the  same. 

When  an  axis  A  passes  through  a  tuple  f^.  then  F  is  a  2/ -palindrome.  Similarly,  if  A  bisects  the  angle 
f 1  then  F  is  a  (2/  +  l)-palindrome. 
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1J3MMA  For  any  k  (0  ^  k  <  m) :  F  is  a  k  -palindrome  iff .k  is  even  and  there  is  an  axis  A  passing 

through  tuple  ffc ,  or  k  is  odd  and  there  is  an  axis  A  bisecting  the  angle  f^ . 


There  are  two  eases  to  consider, 

•  k  is  odd.  Let  k  =  2/  + 1 .  That  F  is  a  (2/  +  l)-palindrome  is  precisely  the  condition  brought  about 
by  A  bisecting  angle  . 

•  k  is  even.  Let  k  =  2/ .  That  F  is  a  2/ -palindrome  is  precisely  the  condition  brought  about  by  A 
passing  through  tuple  . 


We  have  shown  by  the  two  lemmas  that  for  any  k  (0  <  k  <  m),  there  is  an  axis  A  passing  through  tuple  fk  (k 
is  even),  or  an  axis  A  bisecting  angle  (k  is  odd)  precisely  when  there  is  a  match  for  R  in  F ;  at  index  2k .  It 
remains  to  point  out  that  this  is  sufficient  to  capture  all  axes  because  an  axis  A  must  cut  the  circle  at  two  points 
9  apart  and  in  having  k  range  from  0  to  m-1  we  have  covered  a  complete  semi-circle. 

4.  Complexity 

In  this  section  we  show  that  the  complexity  of  the  algorithm  presented  in  Section  2  matches  that  of  the 
problem  and  so  is  optimal,  Le„  the  algorithm  is  0(n  log  n)  and  the  problem  is  shown  to  be  Q(n  log  n).  To 
show  the  latter  we  use  a  reduction  of  the  set  equivalence  problem  to  that  of  deciding  whether  there  is  or  is  not 
an  axis  of  mirror  symmetry  in  a  planar  point  set.  This  also  shows  that  it  is  no  harder  to  find  all  the  axes  of 
mirror  symmetry  than  it  is  to  find  whether  there  are  any  at  all. 

From  the  description  of  the  algorithm  in  section  2  we  can  see  that  it  is  0(n  log  n),  and  this  derives  from  the 
sorting  operation  of  Step  3.  Steps  1,2,4  and  6  are  clearly  0(n)  while  Step  5  can  be  done  in  0(n)  time  using  a 
string  matching  algorithm  such  as  KMP.  The  worst-case  complexity  of  the  KMP  algorithm  is  not  affected  by 
alphabet  size. 

Given  two  sets  of  n  real  numbers  we  wish  to  decide  whether  or  not  they  are  the  same.  This  problem  is 
known  to  be  Q(n  log  n)  work. 

Call  die  sets  A  and  B  and  then  proceed  as  follows: 

1.  Find  the  lowest  value  LA  in  A  and  the  lowest  value  Lg  in  B. 

1  Create  the  set  A 1  from  A  by  adding  the  element  a  -  LA  to  A /  for  each  a  in  A.  Similarly  create  B/ 
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from  B  using  1.B . 

3.  Manufacture  the  2D  point  set  P  defined  as  {(u ,  a ) :  a  €  A 7)  U  {(b ,  -b ) :  b  €  B7}. 

4.  Find  the  axes  of  mirror  symmetry  of  the  set  P.  If  P  has  an  axis  of  mirror  symmetry  then  A  and  B 
are  the  same,  else  they  are  not. 

The  set  P  defined  above  will  have  at  most  one  axis  of  mirror  symmetry  and  this  if  present  will  correspond  to 
the  ’x-axis’  of  the  2D  system  created. 


5.  Rotational  and  point  symmetries 

In  this  section  we  present  a  version  of  the  algorithm  of  section  2  which  will  find  all  the  rotational 
symmetries  of  point  set  P.  Hence,  as  noted  in  the  introduction,  we  will  also  have  discovered  whether  the  set 
has  a  point-symmetry  because  that  is  identical  to  a  rotational  symmetry  of  v  . 


A  rotational  symmetry  is  an  angle  a  s.L  rotation  of  P  by  a  about  its  centroid  is  an  identity  operation  on  P. 
Furthermore,  if  aQ  is  the  smallest  such  angle  then  the  set  of  angles  of  rotational  symmetry  is  precisely  the  set 
of  integer  multiples  of  aQ .  Hence,  it  is  only  necessary  to  find  o0 . 


Given  the  string  representation  F  of  P,  proceed  by  matching  F  against  F 7  .starting  with  the  first  element  of 
F  at  f2  and  moving  to  the  right  This  avoids  a  meaningless  match  of  F  with  itself  as  prefix  of F 7  and  because 
angle  a0  corresponds  to  the  lowest  index  at  which  F  will  match  F7 .  Stop  when  the  first  match  is  found  (or  a 
failure  occurred  at  the  end  of  F 7 ,  implying  no  rotations).  Let  the  index  of  the  match  in  F 7  be  /  ,  which  must 
be  even  as  before.  The  corresponding  aQ  is  the  sum  of  the  angles  in  F7  lying  to  the  left  of  the  match. 


°0  =  fl  +  f3  +  •“  ^-l  • 


The  complexity  of  this  algorithm  is  the  same  as  that  of  the  mirror  symmetry  algorithm.  0(n  log  n). 
Checking  whether  a0  is  an  integer  divisor  of  v  is  constant  time  work  and  so  0(n  log  n)  applies  to  the  point 
symmetry  decision  too. 

We  now  show  that  deciding  whether  P  is  point-symmetric  is  Q(n  log  n)  and  so  the  algorithm  above  for 
rotational  symmetry  and  point  symmetry  is  optimal.  As  in  section  4,  use  the  problem  of  deciding  whether  the 
two  sets  A  and  B  are  identical.  We  show  that  this  is  no  harder  than  deciding  if  P  is  point-symmetric  by  the 
following  algorithm: 

1.  Find  the  lowest  value  LA  in  A  and  the  lowest  value  LB  in  B. 

2.  Create  the  set  A 7  from  A  by  adding  the  element  a  -  LA  to  A 7  for  each  a  in  A.  Similarly  create  B 7 
from  B  using  Lg . 
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3.  Manufacture  the  2*D  point  set  P  defined  as  {(  a ,  a ) :  a  €  A  U  {( ~b ,  -b ) :  b  €  B *}. 

4.  Decide  whether  P  is  point-symmetric,  if  it  is  then  A  and  11  are  the  same,  else  they  are  not 
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